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Summary 

A relatively simple but accurate potential-flow 
model has been developed for studying the 
symmetric- vortex flow on cones. The model is a mod- 
ified version of the model first developed by Bryson, 
in which discrete vortices and straight-line feeding 
sheets were used to represent the flow field. It dif- 
fers, however, in the zero-force condition used to po- 
sition the vortices and determine their circulation 
strengths. The Bryson model imposed the condition 
that the net force on the feeding sheets and discrete 
vortices must be zero. The proposed model satisfies 
this zero-force condition by having the vortices move 
as free vortices, at a velocity equal to the local cross- 
flow velocity at their centers. When the free- vortex 
assumption is made, a solution is obtained in the 
form of two nonlinear algebraic equations that relate 
the vortex-center coordinates and vortex strengths 
to the cone angle and angle of attack. The vortex- 
center locations calculated using the model are in 
good agreement with experimental values. The cone 
normal forces and center locations are in good agree- 
ment with the vortex cloud method of calculating 
symmetric flow fields. 

Introduction 

The symmetric-vortex flow on cones has been 
studied in several investigations in connection with 
the high-angle-of-attack performance of aircraft and 
missiles. Designers would like to use the lift pro- 
duced by these flows to improve maneuvering ca- 
pability, but they also want to understand and be 
able to predict adverse flow patterns, such as vortex 
asymmetry, that can occur. Vortex asymmetries usu- 
ally occur at moderate to high angles of attack and 
can cause serious stability and control problems. A 
thorough knowledge of the symmetric-vortex flow 
field could lead to a better understanding of why and 
under what conditions these asymmetries occur. 

The previous analytical studies of conical vor- 
tex flow have included Navier-Stokes solutions, ob- 
tained mainly at supersonic speeds (refs. 1 to 3), and 
potential-flow solutions, obtained using free-sheet 
methods (refs. 4 to 7). More recently, a vortex cloud 
method (refs. 8 and 9) has been developed that allows 
vortex flow characteristics to be determined for cones 
as well as general forebody shapes. These methods, 
however, involve fairly complicated time-stepping or 
iteration procedures, and the results available in the 
literature cover only a limited number of cone angles 
and angles of attack. The effects of the different pa- 
rameters governing conical vortex flow have not been 
clearly defined. 


This report describes a relatively simple discrete- 
vortex model for studying the vortex flow on cones. 
The model is based on the Bryson model described in 
references 10 and 11, but uses a different zero- force 
condition than that used by Bryson to locate the 
vortices and determine their circulation strengths. 
Bryson imposed the condition that there be zero net 
force on a system consisting of primary vortices and 
straight-line feeding sheets. In the present model, 
this zero-force condition is satisfied by having the 
vortices move as free vortices, that is, at velocities 
equal to the local flow velocities at their centers. 

The assumptions of a free vortex and conical 
flow enable a solution to be obtained in the form 
of two nonlinear algebraic equations that relate vor- 
tex position, vortex strength, and the parameter 
tan a/tan 6. Together with information about the 
separation point, these relationships can be used to 
obtain a solution for a given cone angle and angle of 
attack. 

The vortex centers calculated using the present 
method are shown to be in good agreement with 
the experimental data obtained at incompressible 
laminar-flow conditions in reference 12. When the 
measured separation points are used, the method 
gives an accurate prediction of the travel of the vortex 
centers as the cone goes through an angle-of-attack 
range. The vortex-center locations and normal forces 
predicted using the method also agree with those pre- 
dicted using the vortex cloud method of references 8 
and 9 for the same circulation strength. 

Symbols 

a X-axis influence coefficient 

(symmetric-vortex alignment) 

an,&u X- and Y-axis influence coeffi- 

cients, respectively, that express 
influence of image of vortex 1 on 
conditions at center of vortex 1 for 
asymmetric- vortex alignment (see 
fig. Bl) 

a 12 , b \2 X- and Y-axis influence coefficients, 

respectively, that express influence 
of vortex 2 and its image on con- 
ditions at vortex 1 for asymmetric- 
vortex alignment (see fig. Bl) 

a 2l i ^21 X - and Y-axis influence coefficients, 

respectively, that express influence 
of vortex 1 and its image on con- 
ditions at vortex 2 for asymmetric- 
vortex alignment (see fig. Bl) 


a 22>&22 X- and F-axis influence coeffi- 

cients, respectively, that express 
influence of image of vortex 2 on 
conditions at center of vortex 2 for 
asymmetric-vortex alignment (see 
fig. Bl) 

b Y- axis influence coefficient 

(symmetric-vortex alignment) 

C k nondimensional vortex strength, 

r 

2 xroUoo sin a 

Ck avail nondimensional vortex strength 

made available for vortex formation 
by boundary-layer separation on 
cone 

C kieq nondimensional vortex strength 

required to meet the zero- force 
condition 

Cjy normal- force coefficient, 

cj y local normal- force coefficient, 

measured in X-F plane at a 

given 2 -station along cone axis, 

Local normal force 
(l/2) P U‘iD 

D diameter of cone base 

dp differential pressure across a feeding 

sheet 

F force exerted on vortex system by 

body at a given cross section 

Fp s. pressure force on feeding sheet 

F n local normal force 

F v Kutta-Joukowski force acting at 

vortex, center 

i unit vector along imaginary axis 

k vorticity reduction factor 

L length of cone 

MV momentum of vortex pair 

r radial distance in crossflow plane, 

y/x 2 + y 2 

r 0 radius of cone cross section 

S base area of cone 

t time 

U X-axis component of velocity, 

nondimensionalized with respect 
to {Too sin a, that results from free- 
stream flow around cone at a given 
cross section 


U e velocity at edge of boundary layer 

at a separation point on cone 

«7oo free-stream velocity 

u — X-axis component of velocity 

at a point in flow field, nondimen- 
sionalized with respect to {Too sin a 

V Y -axis component of velocity, 

nondimensionalized with respect 
to {Too sin a, that results from free- 
stream flow around cone at a given 
cross section 

V k induced velocity at vortex center, 

u k + iv k 

v = F-axis component of velocity 

at a point in flow field, nondimen- 
sionalized with respect to {Too sin a 

w complex potential, <j> + iip 

x X-axis coordinate (see fig. 1) 

x f = x/r 0 

y F-axis coordinate (see fig. 1) 

y' = y/r 0 

z Z- axis coordinate (see fig. 1) 

a cone angle of attack 

T circulation, positive clockwise in 

first (+x, -hy) quadrant of X-F 
plane 

r ava ii circulation made available for 

vortex formation by boundary-layer 
separation taking place on cone 

6 cone semi-apex angle 

C = x + fy, complex coordinate of a 

point in flow field 

Co complex coordinate of origin of 

feeding sheet (see fig. Al) 

9 — tan” 1 ^), angle measured from 

positive X-axis to a point in flow 
field 

p fluid density 

4> velocity potential 

<f>s angular distance to separation point 

on cone at a given cross section, 
measured from negative X-axis (see 

fig- 1) 
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^ stream function 

Subscripts and superscripts: 


The real part of equation (1) gives the velocity 
potential 


k primary vortex index 

(”) complex conjugate 

1 refers to conditions at center of 

vortex 1 (see fig. Bl) 

2 refers to conditions at center of 

vortex 2 (see fig. Bl) 


(p = Uqo si n ot 



x l + y z 


x 


r 

27rr 0 f7oo sin a 


tan 1 


( j/jJ/fc \ 

\x-x k ) 


Mathematical Development 

Basic Flow Model 

The model used (see fig. 1) is the two-dimensional, 
discrete-vortex model that is described in refer- 
ence 13 and used by a number of investigators, but 
with a source term added to simulate conical flow. 
The model assumes incompressible potential flow. 
The flow at each cross section is represented by two 
primary vortices placed at the positions x ^ and ±2/^ 
with image vortices placed inside the cone at the po- 
sitions shown in figure 1 to satisfy the tangent-flow 
boundary condition on the surface. The added source 
term has a strength that varies with the distance z 
along the cone axis and simulates the velocities that 
occur normal to the surface at each cross section in 
conical flow. The complex potential for this system 
is 



— r 0 Uoo cos a tan <5 In \J x 2 + y 2 (2) 


The imaginary part of equation (1) gives the 
stream function 


W = <p + iip = Uoo 




In (C - C k) ~ ln 



ip = Uqo si n a 



x z + y l 


y 



In sj(x- X k ) 2 + {y- Vk ) 2 


-i„(c-< k ) + i"(c -g) 

— t 0 Uoo cos a tan 6 ln £ (1) 

The first term on the right-hand side in equa- 
tion (1) is the complex potential for the free-stream 
flow around the cross section; the second through 
fifth terms are those for the primary vortices and 
their images; the last term is the source term added 
for conical flow. This source term assumes the flow 
is purely conical, that is, all flow quantities are con- 
stant along rays emanating from the cone apex. This 
term does not account for the end effects that result 
from the cone having a finite length. (For examples 
of more complex representations that do account for 
end effects, see ref. 14.) 



- In \f(x - x fc ) 2 + (y + Vkf 



— VoUoo cos a tan 6 tan 6 (3) 


The nondimensional velocity in the x direction, 
written in terms of x' and y' coordinates is given by 
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u = r 


1 - 


+ 


V 


x' + y' (x ,2 + y' 2 ) 2 


+ 


27rr 0 Uoo sin a 


W - rfk) 


(x 1 - xjj) 2 + <y' - s')2 


(V-^) 


( x ' - $4) 2 + 0 

{y'+y'k) 

MO 2 

[x'k - X 
1 

f + iy'k + yf 






This velocity can be written as 

tan <5 


U — -f- 


COS0 


tan 


a \j x' 2 + y' 2 


( 5 ) 


where a is an X-axis influence coefficient equal to 

the sum of the terms multiplied by 0 — T }' — , — in 

* Zirr 0 Uoc sm a 

equation (4). 

The velocity in the ^-direction is given by 


This velocity can be written as 


v = V + bCfr + 


tan 6 sin 6 


tan 


a y/ x* 1 + y f 2 


( 7 ) 


where 6 is a F-axis-oriented influence coefficient 
equal^ to the sum of the terms multiplied by 


^roUoosma in equation (6). 


Solution for Conical Flow 

The basic flow problem can be solved if the lo- 
cations and circulation strengths of the primary vor- 
tices are known. Equation (2) defines a velocity po- 
tential <fi that is a function of x and y and satisfies 
the Laplace equation and the boundary conditions 
of the problem. The Laplace equation is automat- 
ically satisfied because each of the singularities in- 
volved satisfies the Laplace equation. On e bound- 
ary condition — that the flow must be tangent to the 
body at the surface — is satisfied by the image system 
used. The other boundary condition — that the flow 
velocity at infinity must be equal to the free-stream 
velocity — is satisfied because the velocity induced by 
each of the singularities varies inversely with the dis- 
tance from the singularity. 

Once this velocity potential is defined, the veloci- 
ties at different points in the potential flow field and 
the forces on the body can be determined. The prob- 
lem then becomes one of determining the locations 
and circulation strengths of the vortex centers for 
different cone angles and angles of attack. These lo- 
cations and strengths are usually determined by im- 
posing a zero-force condition to constrain the motion 
of the primary vortices. 


r 

-( x '~ x 'k) 


2'Kr o U 0Q sin a 

(x'-x' k ) 2 + (y 1 -y 

o 2 


+ (‘'-i s 0 2 


Zero-Force Condition 


The zero-force condition imposed in the present 
method is that the primary vortices move as free 
vortices, that is, at a velocity equal to the local 
fluid velocity of the vortex centers. This condition 
is expressed mathematically by 


v: = 


dtk 

dt 


(8) 


(*'-■'*) V 


X ' k ) 



( x i _ d x > 

r t * 

+ 

CM 

('■ 


2 


+ tan 6 sin0 


This condition is a departure from the zero-force 
condition used by Bryson and others (e.g., see refs. 7 
and 10), for which the force on a feeding sheet that 
extends into the vortex is balanced out by an equal 
and opposite force at the vortex center; thus, there 
is no net force on the vortex system. However, the 
use of equation (8) is justified in appendix A. 
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Conical Flow Assumptions 

Equation (8) relates the induced velocity at the 

vortex center to the absolute velocity ^ of the cen- 
ter relative to the axis-system origin. This absolute 
velocity can be expressed in terms of cone semi- apex 
angle and angle of attack by assuming conical flow 
and making use of the impulsive-flow analogy. Ac- 
cording to the impulsive-flow analogy (e.g., ref. 15), 
the flow at a given cross section of the cone is equiva- 
lent to that about a two-dimensional cylinder whose 
radius is expanding with time, where time in the two- 
dimensional plane is related to the distance along the 
axis of the three-dimensional cone by 


u rr V / 

Uoo cos a 

If the flow is assumed to be conical, the movement of 
the vortex center is known. The center moves radially 
away from the cone axis as either time (in the two- 
dimensional plane) or distance z (along the three- 
dimensional body) increases. Therefore, the velocity 
of the vortex center at a given cross section must 
point in the radial direction and have a magnitude of 

ir = c, “ 00,0 (£) (10) 


Now the derivative with respect to 2 : can be 
written as 

** = (*k) (±i\ = ffe') tan« (11) 
dz \dr 0 J \ dz ) \ dr 0 ) 

Since, for conical flow, the nondimensional distance 
to the vortex center remains constant for a given cone 
angle and angle of attack, 

^ tan 6 Uqo cos a (12) 

dt \r 0 J 

The velocity of the vortex center, nondimensionalized 
with respect to Uoo sin can therefore be written as 

1 f ±k\ = (Tk \ tan<5 ( 13 ) 

Uoo sin a \ dt ) \ r Q ) tan a 


Equations for C k and tan a/tan 6 

In terms of nondimensionalized X- and T-axis ve- 
locity components, the zero-force condition (eq. (8)) 
becomes 

„ tan <5 / cos 8l. \ / r fc \ tan <5 a\ 

+ (14) 

n+wi+r (is) 

K K K tana \r k /r 0 / \r 0 J tana 


(For an extension of the zero-force equations 
(eqs. (14) and (15)) to asymmetric-vortex align- 
ments, see appendix B.) Here, Cfc req is a required 
nondimensional vortex strength — it is needed to meet 
the zero-force condition for a given vortex-center 
location. The influence coefficients a* and bk in 



(see eq. (6)) since, in potential flow calculations, it 
can be assumed that a discrete straight-line vortex 
has no influence on itself. (See ref. 16.) 


Combining the terms containing tan a/tan 6 and 
dividing equation (15) by equation (14) gives 


Yk + h£k 

Uk + G/tCfc 


= tan 9 £ 


(16) 


from which the circulation required to meet the zero- 
force condition can be determined as 

_ -(E/jfctanflfc - Vfc) / 1T) 

*’ req_ (a k twd k -b k ) { 1 


Since the factors defining the nondimensional vortex 
strength in equation (17) — l/fc, V^a*, and b k — are 
functions of x k and y k only, this nondimensional 
strength is uniquely defined by the vortex-center 
location. 

The value of tana/tan<5 corresponding to a given 
vortex-center location and nondimensional vortex 
strength can be determined from either equation (14) 
or (15) as follows: 


tana_ ( r fc/ r o)-(i^fe) cosg fc 
tan 6 (Uk + a kC *) 

tVk + l>„C k ) 

Equations (17) and (18) form the basis of the 
present method; they provide two equations in terms 
of the known parameter tan aj tan 6 and the three un- 
known parameters x r k ,y k , and C k req . A third equa- 
tion that is available for solving the problem is an 
equation for the vortex strength. The nondimen- 
sional vortex strength given by equation (17) is that 
required to meet the zero- force condition, but just 
what value this strength achieves for a given value of 
tan a/tan <5 depends on how much circulation is made 
available for vortex formation by the boundary-layer 
separation that is taking place on the cone. This 
available circulation can be determined from the rate 
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at which circulation is generated at the separation 
points on the cone. 

Determination of Available Circulation 

If the location of the separation point on the cone 
is known, having been either measured or calculated 
using a boundary-layer-separation criterion such as 
that employed in the vortex cloud program of refer- 
ences 8 and 9, then the rate at which circulation is 
generated can be approximated by the relationship: 


avail 

dt 



(19) 


where U e is the edge velocity at the separation point 
(it includes the velocities induced by the vortices). 
The k factor is a vorticity modification factor that 
accounts for the fact that, in general, not all the 
vorticity generated by the boundary- layer separation 
is entrained into the primary vortices. A k factor of 
0.6 is usually assumed. (See ref. 8.) 

Since, for conical flow, T varies linearly with £, 

avail ^avaii 

dt t [ ) 


By using equation (9), equation (19) can be written 
as 


n _ k ( U e \ 2 tana 

avail “ A ( TT . ) T 7 (21) 

47 r Vf/ocsmo/ tan 8 

which is the nondimensional vortex strength available 
for a given value of tan a/tan 6. 


Calculation of Normal-Force Coefficients 

Once the nondimensional vortex strengths and 
locations of the vortex centers have been determined, 
the normal forces on the cone can be calculated from 
the equations derived in appendix C. The equation 
used for calculating the local normal-force coefficient 
at a given cross section is 


C N — 27T 4- 4 — — j— ^ C k sin^J ~ tan<5 sin a cos a (22) 

The equation for the normal-force coefficient for the 
complete cone is 


Cn = 2 




sin dfc 


sin a cos a 


(23) 


Results and Discussion 


Vortex-Center Locations as a Function of 

Cfc and tana/ tan <5 

In the present method, the primary vortices are 
kept force-free by having them move as free vortices, 


that is, at a velocity equal to the local crossflow ve- 
locity at their centers. The vortex-center positions 
and circulation strengths needed to meet this condi- 
tion for different cone angles and angles of attack are 
shown in figure 2. Figure 2 shows contours, obtained 
by iteration from equations (17) and (18), of con- 
stant circulation strength (solid lines) and the values 
of tan a/tan <5 (dashed lines); these contours are plot- 
ted as functions of the coordinates of the primary 
vortex center, x f k and y k . Figure 2 can be used in 
two ways. If the available circulation strength has 
been determined for a given cone angle and angle of 
attack — either by a boundary-layer calculation or by 
experiment — the contours in figure 2 can be used to 
locate the coordinates of the vortex center. If, on 
the other hand, the vortex-center location has been 
measured experimentally, these contours can be used 
to approximate the circulation strengths of the vor- 
tices. In either case, the complete potential-flow field 
can be mapped out, and forces on the body can be 
determined. 

The curves in figure 2 show a limit to vortex- 
center travel. The vortex centers would not be ex- 
pected to go beyond the curve for which tan a /tan 6 
is infinite (a cone angle of attack of 90°). The curve 
that defines this limit is the so-called Foppl curve. 
(See refs. 13 and 17.) The Foppl curve gives the 
locations of the vortex centers for which the veloc- 
ity V£ at the center is zero for a constant-diameter 
cylinder in two-dimensional flow. It can be shown 
that this condition corresponds to having the denom- 
inators in the equations for tan a/tan 6 (eqs. (18)) 
equal to zero, which would produce an infinite value 
of tan a/tan 5. The Foppl curve, therefore, defines 
the limit of vortex-center travel for conical flow. 

Comparison of Calculated Center 
Locations With Experiment 

Figure 3 shows measured vortex-center locations 
from reference 12 superimposed on the constant C * 
and tan a/tan 8 curves of figure 2. The experimental 
results were obtained in a water tunnel at a Reynolds 
number (based on cone length) of about 2.7 x 10 4 . 
(The boundary layers for this condition were de- 
scribed in ref. 12 as being steady and laminar.) The 
results are shown for cone semi-apex angles of 7.5° 
and 12.5° in figures 3(a) and 3(b), respectively. The 
values of tan a/tan 6 for individual data points are 
indicated in the figures. These results show that as 
the jingle of attack is increased, the vortex center 
moves upward, toward curves for higher values of 
Q., and to the right, toward curves for higher val- 
ues of tan a/tan 5. For both cones, the measured 
vortex-center location for a given value of tan a/tan 6 
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is approximately that needed to meet the zero-force 
requirement of the present method. At a value of 
tan a/tan 8 of approximately 2.0, for example, the 
experimental data indicate a vortex-center location 
at about x^jr 0 = 1.16 and ykjr 0 = 0.33; these val- 
ues are in good agreement with those needed to meet 
the zero- force requirement. 

The nondimensional vortex strengths of the pri- 
mary vortices for the test results shown in figure 3 
were not defined in reference 12; however, they were 
calculated in the present investigation by using the 
measured separation point locations given in refer- 
ence 12. These measured locations are shown for 
the 7.5° and 12.5° cones in figure 4(a). The vortex 
strengths for a given value of tan a / tan 6 and a given 
separation point were determined by constructing, 
for that value of tan a/tan <5, a curve of yk/r Q ver- 
sus Xk/r 0 similar to those shown in figures 2 and 
3, and then iterating along this curve until the vor- 
tex strength required to meet the zero-force condition 
(eq. (17)) matched that being made available at the 
separation point (eq. (21)). The nondimensional vor- 
tex strengths determined by this iteration procedure 
for the two cones are presented in figure 4(b), which 
shows that this vortex strength was approximately a 
linear function of tan a/tan 6. 

Figures 5(a) and 5(b) compare the nondimen- 
sional vortex strengths given in figure 4(b) with those 
determined by using the Bryson equations; these fig- 
ures also compare the vortex-center locations deter- 
mined by each method with the experimental-center 
locations obtained from reference 12. The vortex- 
center locations calculated using the present method 
are in generally good agreement with the measured 
values over the range of values of tan a/tan 6 shown. 
The predictions made with the Bryson method, on 
the other hand, were not nearly as good. Because 
of the relationship between (fr s and tan a/tan 6 con- 
tained in the Bryson zero-force equation, solutions 
with the Bryson method were possible for only cer- 
tain segments of the curves of (f) s versus tan a/tan 6 
shown in figure 4(a). That is, for the 7.5° cone, solu- 
tions were possible only for values of tan a/tan <5 be- 
tween about 3.6 and 4.4; for the 12.5° cone, no solu- 
tions were possible over the complete range of values 
of tana/tan 8 for which data were taken. The nondi- 
mensional vortex strengths determined by the Bryson 
method were significantly higher than those deter- 
mined by the present method. The vortex-center 
locations were much farther away from the axis of 
asymmetry for a given value of tan a/tan 6 than in- 
dicated by either the present calculations or the ex- 
perimental data. 


Comparison With Vortex Cloud Results 

To further verify the basic assumption of the 
present method — for the cone, the vortices remain 
force-free by traveling with the local crossflow 
velocity — calculations made using the present 
method are compared in this section with predictions 
made using the vortex cloud method of references 8 
and 9. The vortex cloud method uses discrete vor- 
tices generated at the separation points on the body 
to model the vortex wake. These discrete vortices are 
assumed to travel at the local crossflow velocity and 
eventually wrap up into the equivalent of the primary 
vortices modeled in the present investigation. The 
method uses a modified Stratford separation criterion 
to determine the locations of the separation points 
and assumes circulation is generated at the rate in- 
dicated by equation (19) at these points. To provide 
a common basis of comparison for the two methods, 
the calculations for both are based on the separation- 
point locations and nondimensional vortex strengths 
determined by the vortex cloud method. 

The flow fields depicted by the two methods are 
compared in figure 6. The points indicated by the 
symbols are the locations of the individual vortex el- 
ements used in the vortex cloud method. The stream- 
lines indicated in figure 6 are those calculated using 
the present method. (The streamlines shown are not 
actually streamlines, but rather, projections in the 
crossflow plane of the paths followed by individual 
particles in the three-dimensional flow.) The vortex- 
center locations indicated by the two methods are in 
good agreement for the three angles of attack shown. 

Normal-force coefficients calculated using the 
present method and the vortex cloud method are 
compared in figures 7 to 9. Figure 7 shows the 
variations of local normal- force coefficient c ^ with 
distance along the body z/L for different angles of 
attack. Figures 8 and 9 show the variations in to- 
tal normal- force coefficient with tan a/tan 8 for cone 
semi-apex angles of 5° and 7.5°, respectively. The 
nondimensional vortex strengths used in obtaining 
these results were those determined for the different 
cone angles and angles of attack by the vortex cloud 
program; these strengths are shown as a function of 
tan a/tan 6 in figures 8(a) and 9(a). 

Figure 7 shows generally good agreement between 
the local normal-force coefficients calculated using 
the two methods, based on the same nondimensional 
vortex strength. The local normal-force coefficients 
calculated using the vortex cloud method show some 
oscillation with increasing values of z/L, particularly 
for the 5° cone. This oscillation is probably the result 
of errors in the iteration procedure used in the vortex 
cloud calculations for these particular cases. The 
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present method shows a straight-line variation in 
with z/L , which is in keeping with the conical-flow 
assumptions on which the method is based. 

To indicate the level of the normal force produced 
by vortex flow on the two cones, figures 8 and 9 
show the attached-flow normal-force coefficients (the 
normal- force coefficients for the cones without vor- 
tex flow, as determined from slender-body theory) 
and the total normal-force coefficients (the sum of 
the attached-flow and vortex-flow normal-force coef- 
ficients). Figures 8 and 9 show that the nondimen- 
sional vortex strengths — as calculated by the vortex 
cloud method— were slightly higher for the 7.5° cone 
than for the 5° cone at a given value of tan a/tan 6. 
These higher values of resulted in a noticeably 
higher level of total normal force for the 7.5° cone. 
The normal-force coefficients for the vortex cloud 
method were generally lower than those calculated 
by the present method. At least some of this lower 
level of overall normal force was the result of a lower 
level of local normal force produced toward the end 
of the cone by end effects. The overall level of agree- 
ment, however, was good. 

Conclusions 

A potential-flow model has been developed for 
studying the symmetric-vortex flow on cones. The 
model assumes conical flow conditions and satisfies 
the zero-force condition on the primary vortices by 


having them move at the local induced-flow veloc- 
ity. The solution obtained is in the form of two 
nonlinear algebraic equations that relate the vortex- 
center locations, the vortex strengths, and the ratio 
of the tangent of the angle of attack to the tangent of 
the cone semi-apex angle tana/tan 6. The following 
conclusions can be derived from the results of this 
investigation: 

1. The assumption that the vortex maintains a zero- 
force condition by being a free vortex, that is, by 
moving at a velocity equal to the local crossflow 
velocity, appears valid. 

2. When measured separation-location points were 
used to calculate the available circulation, the 
vortex-center locations calculated using the 
present method were in good agreement with low- 
speed laminar water-tunnel data. 

3. The vortex-center locations and normal forces cal- 
culated using this method were in good agree- 
ment with those calculated using the vortex cloud 
method of references 8 and 9 for the same vortex 
strengths. 

NASA Langley Research Center 
Hampton, VA 23665-5225 
March 22, 1990 
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Appendix A 

Justification for Zero-Force Condition 
Used in Present Method 

Zero-Force Condition Used by Bryson 

The zero- force condition used to determine vor- 
tex positions and strengths in a concentrated vortex 
model is based on the principle that an isolated po- 
tential vortex should sustain no force. (See ref. 7.) In 
the Bryson method (ref. 10), the zero- force equation 
is developed by installing a feeding sheet (fig. Al) and 
balancing out the force on the feeding sheet with an 
equal and opposite force on the vortex. This feeding 
sheet is assumed to be the path along which vor- 
ticity is fed into the primary vortex and acts as a 
mathematical branch cut, which is needed to make 
the velocity potential for the system single-valued. 
There is a discontinuity in <j> equal to the circulation 
strength across this branch cut, and since T varies 
with distance along the cone axis, </> also varies with 
this distance. Since time and distance along the cone 
axis are related (see eq. (9) of text), there will be, ac- 
cording to the unsteady Bernoulli equation, a pres- 
sure difference across the feeding sheet of 

dp = p ^ (Al) 

When integrated over the length of the feeding sheet, 
this pressure difference results in a pressure force of 

f f.s. = (Cfc - Co) (A2) 

Bryson required that this force be balanced out by 
an equal and opposite force at the vortex center. 

The force at the vortex center arises from a non- 
zero local-flow velocity at the center, and according 
to the Kutta-Joukowski theorem, is given by 

F '= ! ' ,r ( v i'-f) < as > 

This force consists of two parts: (1) the force ipTV£ 
is due to the induced effects at the vortex center. 
(Here, V£, the induced velocity, is equal to the free- 
stream crossflow velocity plus the velocities induced 
by all of the singularities in the system, except the 
primary vortex itself; the velocity induced by a pri- 
mary vortex at its own center is assumed to be zero.) 

(2) The force —ipT^ is due to the absolute motion 
of the vortex center. 


Combining equations (A2) and (A3) gives the 
following zero- force equation used by Bryson: 

-vf - <■>> + it,T ( v * ~ w) = 0 (A4) 

Zero-Force Condition Used in Present 

Method 

The zero-force equation used by Bryson repre- 
sents one way of resolving the problem of how to 
make the velocity potential single- valued and at the 
same time ensure zero net force on the vortex system. 
This solution is not completely satisfactory, however, 
since it requires a pressure force on the feeding sheet 
that is balanced out only “in the mean” to get zero 
net force on the vortex system. In the present paper, 
a new approach is taken; the branch cut is assumed 
to lie between the primary vortex and its image (in- 
stead of between the primary vortex and a separa- 
tion point). Then, when the equation for the forces 
on the vortex system is written out, it shows that 
the pressure force on the branch cut is actually part 
of the change in momentum that the vortex system 
undergoes in response to the force exerted on it by 
the body. Therefore, the pressure force should not 
be balanced out. (See fig. (A2).) 

The momentum of a pair of vortices of equal cir- 
culation strength is equal to the product of the fluid 
density, the circulation, and the distance between the 
vortices. (See refs. 18 and 19.) For the system shown 
in figure (A2), the total momentum is given by the 
following: 

(A5 > 

The force F that the body exerts on the fluid is 
equal to the rate of change of momentum. Taking the 
derivative of MV with respect to time and expanding 
the results gives 

(I) <A«> 

If the branch cuts are taken between the primary 
vortices and their images, as shown in figure (A2), 
the first term in equation (A6) can be interpreted 
as a force on the upper branch cut and the sec- 
ond term as a force on the lower branch cut. These 
forces are of the same form as Bryson’s feeding-sheet 


MV 


= ip T ^Ofc 
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force (eq. (A2)), but they act over the distance be- 
tween the primary and image vortices instead of the 


length - £ 0 . The terms ipT^ and -ipF~j^ in 
equation (A6) can be interpreted as forces caused 
by the absolute motion of the primary vortices; like- 
wise. the terms — f 


ipT & (a) and ipT ® (§) 


can 


be interpreted as forces associated with the motion 
of the image vortices. (These forces are shown as 
dashed vectors in fig. (A2) because, in accordance 
with d’Alembert’s principle, they are effective forces, 
equal but opposite in direction to rates of change in 
momentum, rather than actual forces.) 

The effective forces shown in figure (A2) combine 
with the force F that the body exerts on the fluid to 
form a system in equilibrium. Since these forces are 
in equilibrium, there cannot be an additional force 
that is due to induced effects acting at the vortex 
center; if there were, the force on the body would not 
be equal to that given by the momentum theorem. 

The zero-force condition to impose with the vor- 
tex system in figure (A2), therefore, is that there can 
be no force due to induced effects acting at the vortex 
center. For the cone in the present study, this con- 
dition was imposed by requiring that each primary 
vortex move as a free vortex, that is, at a velocity 
equal to the local flow velocity at its center: 


vi = 


d tk 

dt 


(A7) 


In the present method, the feeding sheet has es- 


sentially been removed from consideration; a feeding 
sheet per se is not considered necessary with this type 
of concentrated vortex model. The branch cut from 
the primary vortex center to its image vortex center 
serves to make the velocity potential single- valued. 
And while the vortices are assumed to be fed in some 
way so that their circulation strengths increase with 
distance along the cone axis, the actual manner in 
which they are fed is not considered important. In 
the actual flow, vorticity is fed into primary vortices 
along curved feeding sheets, as best depicted in the 
vortex cloud or the free vortex sheet methods. In 
those methods, however, as in the present method, 
there are no pressure forces associated with the feed- 
ing sheets that would have to be balanced out by 
a force on the vortex. The curved sheets in thosb 
methods are assumed to be free surfaces that cannot 
sustain a pressure force. 

Also, the requirement that there be no force due 
to induced effects at the vortex center does not 
necessarily mean that the vortex moves with the local 
crossflow velocity in all types of flow fields. In an 
unsteady flow field, the vortex velocity must be equal 
to the local flow velocity at a given instant; however, 
since the vortex can also undergo motions caused by 
unsteady effects, the vortex velocity, in general, is not 
equal to the local flow velocity. The requirement of 
equation (A7), therefore, does not necessarily apply 
to problems such as the vortex flow over cylinders 
and narrow delta wings, where other constraints on 
the vortex system may have to be considered. 
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Appendix B 

Extension of Zero- Force Condition to 
Asymmetric Vortices 

The zero- force condition on which the present 
method is based — that the primary vortices move 
as free vortices— can be extended to include vor- 
tices that have different circulation strengths and are 
asymmetrically aligned about the y = 0 axis. (See 
fig. Bl.) This is done by writing a zero- force equa- 
tion similar to equations (14) and (15) at each of the 
two primary vortex centers as follows: 

At vortex 1, 


U\ 4 anCfci 4 a\2C k 2 = 

(s- 

1 

ri/r 0> 

\ tan 8 

I cos 8 1 

/ tana 

(Bl) 

V\ 4 6 n C 7 *x 4 b\2C k 2 = | 

(s- 

1 ) 
r i/r„ J 

i tan 6 . 

sin 81 

' tana 

(B2) 

At vortex 2 , 





U2 4 o.2\Ck\ 4 022^2 = 

(5- 

2 > 
r 2 / T o / 

\ tan 6 

) cos 82 

! tana 

(B3) 

V2 4 b2\Ck\ 4 b22Ck2 = | 

( ^2 _ 
\ T ° 

1 ) 
T 2 /r 0 ) 

tan 6 , 

sin 82 

tan a 

(B4) 


Combining equations (Bl) and (B2) gives the 
following equation: 

(an tan 0i — b\\)C k \ + (ai2 tan0i 

- &12) Qb2 = - {U\ tan 0i - Vi) (B5) 


Combining equations (B3) and (B4) gives the 
following equation: 

(a2i tan 02 - b 2 i)C k i + (a 2 2 tan 02 

— 622) 0*2 = - (U 2 tan 02 - V 2 ) (B6) 

Equations (B5) and (B6) can be solved simultane- 
ously to give the circulation strengths Cl and C2 for 
assumed primary vortex positions. Equations (Bl) 
and (B3) can then be solved to give the following val- 
ues of tan aj tan 6 (designated by subscripts 1 and 2 
for vortices 1 and 2, respectively). 

From equation (Bl), 

f = ( ri/r °~ nfe ) cos<?1 (B7) 

\tan 6/1 U\ + a\iC k i + ai 2 C k 2 

From equation (B3), 

/tana\ = {r 2 /r 0 - cosfl 2 

\tan6/ 2 U 2 + d 2 \C k \ + & 22 C k 2 ' 

A solution is obtained when the vortex centers are 
located such that 

/ tan a\ ( tan a \ , 

\ tan 6 / j = ( tan 6/2 ^ 

Asymmetric solutions have been obtained using 
the basic Bryson equations in reference 20. Solutions 
using the general asymmetric form of the present 
method have been obtained in reference 21. 
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Figure Bl. Asymmetric- vortex alignment at a cross section. 
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Appendix C 

Derivation of Equations for Cone Normal 
Force 

The normal force for the cone with vortex flow is 
determined from the rate of change in the momentum 
of the vortex system at a given cross section. The 
momentum of a pair of vortices is equal to the fluid 
density times the product of their circulation and 
the distance between them. (See refs. 18 and 19.) 
Since two pairs of vortices are involved — two primary 
vortices and their images — the total momentum for 
the vortex system at a given cross section can be 
written, in terms of the real variables r k and r 0 , as 


MV = 2 pT r fc - — 
\ r k 


where the momentum is directed normal to the line 
between the image and the primary vortex. The 
normal force is the rate of change of this momentum 
in the x-direction and is given by 

F " = 2p {i)l r ( rt ~n) sin5t (C2) 

Taking the derivative with respect to t gives 


F n — 2p 


) ( l - 

r 0 r k /r 0 J \ 


r o r r o . 

it T+ -dr) mne ’’ 


Since r Q ,r k , and T vary linearly with time, the fol- 
lowing relationships can be substituted into equa- 
tion (C3): 


dr Q dr a ( dz\ 

~df = ~dl{di) :=Uoctan6cOSa (C5) 


dr r r . . 

= — = — (Uoo) tan 6 cos a ( 

The local normal- force coefficient then becomes 

_ Fn 167Tf 0 f 

(1/2 )pUlD D Vr-o 


C k tan 6 sin a cos a sin 0 k (C7) 


r 0 ztax\6 z 
D ~ 2ltan5 = 2 L 

equation (C7) can be written as 
(r k l \ z 

C N — o7r f — j -~C k tan 6 sin a cos a sin 9 k 

\ r o r k! r oJ 

(C9) 

When equation (C9) is added to the usual (slender- 
body) value of local normal-force coefficient, the 
equation for local normal-force coefficient becomes 


Cjv = 2n 1 -f 4 


— jCjtsin^ j tan 6 sin a cos a (CIO) 


The normal-force coefficient for the complete cone is 


Cv = — J c N dz = 2 j^l 

+ 4 - — y— ^ Ck sin#*. sin a cos a (Cl 1) 

which is equivalent to the equation used by Bryson 
in reference 10. 
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(b) Vortex-system geometry in crossflow plane. 
Figure 1 . Schematic of vortex-flow model. 











O Experiment (ref. 11) 

Present method 

Bryson method (ref. 20) 




(a) 5 = 7.5°. (b) 5 = 12.5°. 

Figure 5. Nondimensional vortex strengths and vortex-center locations determined from experimental separa- 
tion points given in figure 4(a). 
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z/L 

(a) 6 = 5°. 

Figure 7. Comparison of local normal-force coefficients calculated using present method and vortex cloud 
method. 
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(b) <5 = 7.5°. 
Figure 7. Concluded. 





tan a/tan 5 


(b) Cjy versus tan a/ tan 6. 

Figure 9. Comparison of total normal-force coefficients calculated using present method and vortex cloud 
method. <5 = 7.5°. 




NASA 

National Aeronautics and 
Space Administration 


1. Report No. 

NASA TP-2989 

4. Title and Subtitle 


Report Documentation Page 


2. Government Accession No. 


Discrete- Vortex Model for the Symmetric- Vortex Flow on Cones 


7. Author(s) 

Thomas G. Gainer 


9. Performing Organization Name and Address 

NASA Langley Research Center 
Hampton, VA 23665-5225 


3. Recipient’s Catalog No. 


5. Report Date 

May 1990 


6. Performing Organization Code 


8. Performing Organization Report No. 

L- 16586 


10. Work Unit No. 

505-60-21-06 

11. Contract or Grant No. 


12. Sponsoring Agency Name and Address 

National Aeronautics and Space Administration 
Washington, DC 20546-0001 


13. Type of Report and Period Covered 

Technical Paper 


14. Sponsoring Agency Code 


15. Supplementary Notes 


16. Abstract 


A relatively simple but accurate potential-flow model has been developed for studying the 
symmetric- vortex flow on cones. The model is a modified version of the model first developed by 
Bryson, in which discrete vortices and straight-line feeding sheets were used to represent the flow 
field. It differs, however, in the zero-force condition used to position the vortices and determine 
their circulation strengths. The Bryson model imposed the condition that the net force on the 
feeding sheets and discrete vortices must be zero. The proposed model satisfies this zero-force 
condition by having the vortices move as free vortices, at a velocity equal to the local crossflow 
velocity at their centers. When the free-vortex assumption is made, a solution is obtained in the 
form of two nonlinear algebraic equations that relate the vortex-center coordinates and vortex 
strengths to the cone angle and angle of attack. The vortex-center locations calculated using 
the model are in good agreement with experimental values. The cone normal forces and center 
locations are in good agreement with the vortex cloud method of calculating symmetric flow fields. 


17. Key Words (Suggested by Authors(s)) 
Aeronautics 
Aerodynamics 


18. Distribution Statement 

Unclassified — Unlimited 


Subject Category 02 


19. Security Classif. (of this report) 

Unclassified 


20. Security Classif. (of this page) 

Unclassified 


NASA FORM 1626 oct 86 

For sale by the National Technical Information Service, Springfield, Virginia 22161-2171 


NAS A- Langley, 1990 

















iiiiii mm 





